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Abstract —This letter proposes a parametric sparse multi¬ 
ple input multiple output (MIMO)-OFDM channel estimation 
scheme based on the finite rate of innovation (FRI) theory, 
whereby super-resolution estimates of path delays with arbitrary 
values can be achieved. Meanwhile, both the spatial and tem¬ 
poral correlations of wireless MIMO channels are exploited to 
improve the accuracy of the channel estimation. For outdoor 
communication scenarios, where wireless channels are sparse in 
nature, path delays of different transmit-receive antenna pairs 
share a common sparse pattern due to the spatial correlation 
of MIMO channels. Meanwhile, the channel sparse pattern is 
nearly unchanged during several adjacent OFDM symbols due to 
the temporal correlation of MIMO channels. By simultaneously 
exploiting those MIMO channel characteristics, the proposed 
scheme performs better than existing state-of-the-art schemes. 
Furthermore, by joint processing of signals associated with 
different antennas, the pilot overhead can be reduced under the 
framework of the FRI theory. 

Index Terms —Super-resolution, sparse channel estimation, 
MIMO-OFDM, finite rate of Innovation (FRI), spatial and 
temporal correlations. 

I. Introduction 

M ultiple input multiple output (MIMO)-OFDM is 
widely recognized as a key technology for future 
wireless communications due to its high spectral efficiency 
and superior robustness to multipath fading channels fll . 

For MIMO-OFDM systems, accurate channel estimation is 
essential to guarantee the system performance El. Gener¬ 
ally, there are two categories of channel estimation scheme 
for MIMO-OFDM systems. The first one is nonparametric 
scheme, which adopts orthogonal frequency-domain pilots 
or orthogonal time-domain training sequences to convert the 
channel estimation in MIMO systems to that in single antenna 
systems El. However, such scheme suffers from high pilot 
overhead when the number of transmit antennas increases. 
The second category is parametric channel estimation scheme. 
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which exploits the sparsity of wireless channels to reduce 
the pilot overhead El^ 01 . The parametric scheme is more 
favorable for future wireless systems as it can achieve higher 
spectral efficiency. However, path delays of sparse channels 
are assumed to be located at the integer times of the sampling 
period, which is usually unrealistic in practice. 

In this letter, a more practical sparse MIMO-OFDM channel 
estimation scheme based on spatial and temporal correlations 
of sparse wireless MIMO channels is proposed to deal with 
arbitrary path delays. The main contributions of this letter are 
summarized as follows. First, the proposed scheme can achieve 
super-resolution estimates of arbitrary path delays, which is 
more suitable for wireless channels in practice. Second, due 
to the small scale of the transmit and receive antenna arrays 
compared to the long signal transmission distance in typical 
MIMO antenna geometry, channel impulse responses (CIRs) 
of different transmit-receive antenna pairs share common path 
delays Q, which can be translated as a common sparse pattern 
of CIRs due to the spatial correlation of MIMO channels. 
Meanwhile, such common sparse pattern is nearly unchanged 
along several adjacent OFDM symbols due to the temporal 
correlation of wireless channels 0, Q. Compared with 
previous work which just simply extends the sparse channel 
estimation scheme in single antenna systems to that in MIMO 
by exploiting the spatial correlation of MIMO channels Q 
or only considers the temporal correlation for single antenna 
systems 0, Q, the proposed scheme exploits both spatial 
and temporal correlations to improve the channel estimation 
accuracy. Third, we reduce the pilot overhead by using the 
finite rate of innovation (FRI) theory 0, which can recovery 
the analog sparse signal with very low sampling rate, as a 
result, the average pilot overhead per antenna only depends 
on the channel sparsity level instead of the channel length. 

Notation: (.)t and {■)^ are the Moore-Penrose matrix 
inversion operation and matrix conjugate transpose operation, 
respectively. diag{x} is a diagonal matrix with the vector x 
on its diagonal. The operator * denotes the linear convolution. 

11. Sparse MIMO Channel Model 

The MIMO channel is shown in Fig. [H and its following 
characteristics will be considered in this letter. 

1) Channel Sparsity: In typical outdoor communication 
scenarios, the CIR is intrinsically sparse due to several sig¬ 
nificant scatterers 0, 0. For an Nt x Nr MIMO system, 
the CIR (f) between the ith transmit antenna and the jth 
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Fig. 1. Spatial and temporal correlations of wireless MIMO channels. 


receive antenna can be modelled as ID, 

= E l<i<Nt,l<j<Nr, (1) 

p^l 

where i5(-) is Dirac function, P is the total number of resolv¬ 
able propagation paths, Tp’^'^ and ap’^'^ denote the path delay 
and path gain of the pth path, respectively. 

2 ) Spatial Correlation: Because the scale of the transmit or 
receive antenna array is very small compared to the long signal 
transmission distance, channels of different transmit-receive 
antenna pairs share very similar scatterers. Meanwhile, for 
most communication systems, the path delay difference from 
the similar scatterer is far less than the system sampling period. 
Therefore, CIRs of different transmit-receive antenna pairs 
share a common sparse pattern, although the corresponding 
path gains may be quite different 0. 

3) Temporal Correlation: For wireless channels, the path 
delays vary much slowly than the path gains, and the path 
gains vary continuously 0. Thus, the channel sparse pattern 
is nearly unchanged during several adjacent OFDM symbols, 
and the path gains are also correlated 0. 


III. Sparse MIMO-OFDM Channel Estimation 

In this section, the widely used pilot pattern is briefly 
introduced at first, based on which a super-resolution sparse 
MIMO-OFDM channel estimation method is then applied. 
Finally, the required number of pilots is discussed under the 
framework of the FRI theory. 

A. Pilot Pattern 

The pilot pattern widely used in common MIMO-OFDM 
systems is illustrated in Fig. |2] In the frequency domain, Np 
pilots are uniformly spaced with the pilot interval D (e.g., 
D = 4 in Fig. |2|l. Meanwhile, every pilot is allocated with a 
pilot index I for 0 < I < Np — 1, which is ascending with the 
increase of the subcarrier index. Furthermore, to distinguish 
MIMO channels associated with different transmit antennas, 
each transmit antenna uses a unique subcarrier index initial 
phase 9i for 1 < i < Nt and {Nt — i)Np zero subcarriers to 
ensure the orthogonality of pilots 0. Therefore, for the ith 
transmit antenna, the subcarrier index of the Ith pilot is 

= + 0</<iVp-i. (2) 

Consequently, the total pilot overhead per transmit antenna 
is iVp_totai = NfNp, and thus Np can be also referred as the 
average pilot overhead per transmit antenna in this letter. 
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Fig. 2. Pilot pattern. Note that the specific Nt = 2, D = 4, Np = 4, 
Pip total = 8 are used for illustration purpose. 


B. Super-Resolution Channel Estimation 

At the receiver, the equivalent baseband channel frequency 
response (CFR) H{f ) can be expressed as 


H{f) = E - /,/2 < / < /«/2, (3) 

where superscripts i and j in o are omitted for conve¬ 
nience, fs = 1/Ts is the system bandwidth, and Ts is the 
sampling period. Meanwhile, the iV-point discrete Fourier 
transform (DFT) of the time-domain equivalent baseband 
channel can be expressed as 0, i.e., 

H[k]=H{^), 0<k<N-l. (4) 


Therefore, for the (i,j)th transmit-receive antenna pair, ac¬ 
cording to (|2|l - (IHl, the estimated CFRs over pilots can be 
written as 
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p=l 


= = Hi 
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(5) 


where for 0 < ^ < Np — 1 can be obtained by using 

the conventional minimum mean square error (MMSE) or least 
square (LS) method Cl, and is the additive white 

Gaussian noise (AWGN). 

Eq. © can be also written in a vector form as 


where 

Si 

g-f27r4f ^ 

Because the wireless channel is inherently sparse and the 
small scale of multiple transmit or receive antennas is neg¬ 
ligible compared to the long signal transmission distance, 
CIRs of different transmit-receive antenna pairs share common 
path delays, which is equivalently translated as a common 
sparse pattern of CIRs due to the spatial correlation of 
MIMO channels 0, i.e., Tp’^^ = Tp and = v[Z] 

for I < p < P, 1 < i < Nt, 1 < j < Np. Hence by exploiting 
such spatially common sparse pattern shared among different 
receive antennas associated with the ith transmit antenna, we 
have 
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H* = VA* H- W\ 1 < i < Nt, (7) 

where the Np x Np measurement matrix H* is 

7jr(i.iVr)[o] 
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V = [v[0], v[l], • • • ,v[Np — 1]]^ is a Vandermonde matrix of 
size Np X Np, A* = • ,a(*’'^’')] of size Np x 

Nr, and W® is an Np x Nr matrix with in its jth 

column and the {I + l)th row. 

When all Nt transmit antennas are considered based on Q, 
we have 


H = VA + W, (8) 

where H = - ,H^*] of size Np x NtNr, A = 

[Ai,A2,--- ,A^‘], and W = ,W^‘]. 

Comparing the formulated problem ® with the classical 
drrection-of-arrival (DOA) problem Q, we find out that they 
are mathematically equivalent. Specifically, the traditional 
DOA problem is to typically estimate the DOAs of the P 
sources from a set of time-domain measurements, which are 
obtained from the Np sensors outputs at NtNr distinct time 
instants (time-domain samples). In contrast to our problem 
in (Ell, we try to estimate the path delays of P multipaths from 
a set of frequency-domain measurements, which are acquired 
from Np pilots of NtNr distinct antenna pairs (antenna- 
domain samples). It has been verified in lITOll that the total least 
square estimating signal parameters via rotational invariance 
techniques (TLS-ESPRIT) algorithm in ||9l can be applied 
to ([8]l to efficiently estimate path delays with arbitrary values. 

By using the TLS-ESPRIT algorithm, we can obtain super¬ 
resolution estimates of path delays, i.e., fp for 1 < p < P, 
and thus V can be obtained accordingly. Then path gains can 
be acquired by the LS method Q, i.e., 

A = VtH = (V^V)-iV^H. (9) 


Eor a certain entry of A, i.e., dpbecause 6i is known 
at the receiver and Tp has been estimated after applying the 
TLS-ESPRIT algorithm, we can easily obtain the estimation 
of the path gain dpfor l<p<P, !<*< Nt, 1 < 
j < Nr- Einally, the complete CER estimation over all OEDM 
subcarriers can be obtained based on (O and (|4|i. 

Eurthermore, we can also exploit the temporal correlation 
of wireless channels to improve the accuracy of the channel 
estimation. Errst, path delays of CIRs during several adjacent 
OEDM symbols are nearly unchanged 0, Q, which is 
equivalently referred as a common sparse pattern of CIRs 
due to the temporal correlation of MIMO channels. Thus, 
the Vandermonde matrix V in (|8]l remains unchanged across 
several adjacent OEDM symbols. Moreover, path gains during 
adjacent OEDM symbols are also correlated owing to the 
temporal continuity of the CIR, so A’s in (l8]l for several 
adjacent OEDM symbols are also correlated. Therefore, when 
estimating CIRs of the qth OEDM symbol, we can jointly 
exploit H’s of several adjacent OEDM symbols based on ([8]l, 
i.e.. 
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( 10 ) 


where the subscript p is used to denote the index of the 
OEDM symbol, and the common sparse pattern of CIRs is 
assumed in 2R + 1 adjacent OEDM symbols ||7]. In this way. 


the effective noise can be reduced, so the improved channel 
estimation accuracy is expected. 

In contrast to the existing nonparametric scheme which 
estimates the channel by interpolating or predicting based 
on CERs over pilots m, El, our proposed scheme exploits 
the sparsity as well as the spatial and temporal correlations 
of wireless MIMO channels to first acquire estimations of 
channel parameters including path delays and gains, and then 
obtain the estimation of CER according to ([^i and (|4]i. 

C. Discussion on Pilot Overhead 

Compared with the model of the multiple filters bank based 
on the ERI theory ifTOll . it can be found out that CIRs of NtNr 
transmit-receive antenna pairs are equivalent to the NtNr 
semi-period sparse subspaces, and the Np pilots are equivalent 
to the Np multichannel filters. Therefore, by using the ERI 
theory, the smallest required number of pilots for each transmit 
antenna is Np = 2P in a noiseless scenario. Eor practical 
channels with the maximum delay spread Tmax, although the 
normalized channel length L = rmax/Ts is usually very large, 
the sparsity level P is small, i.e., P « L 0. Consequently, 
in contrast to the nonparametric channel estimation method 
where the required number of pilots heavily depends on L, our 
proposed parametric scheme only needs 2P pilots in theory. 
Note that the number of pilots in practice is larger than 2P to 
improve the accuracy of the channel estimation due to AWGN. 

IV. Simulation Results 

A simulation study was carried out to compare the per¬ 
formance of the proposed scheme with those of the exist¬ 
ing state-of-the-art methods for MIMO-OEDM systems. The 
conventional comb-type pilot and time-domain training based 
orthogonal pilot (TTOP) 0 schemes were selected as the typi¬ 
cal examples of the nonparametric channel estimation scheme, 
while the recent time-frequency joint (TEJ) channel estimation 
scheme 11 was selected as an example of the conventional 
parametric scheme. System parameters were set as follows: 
the carrier frequency is /c = 1 GHz, the system bandwidth 
is /s = 10 MHz, the size of the OEDM symbol is A = 4096, 
and Ng = 256 is the guard interval length, which can 
combat channels whose maximum delay spread is 25.6 p,s. The 
International Telecommunication Union Vehicular B (ITU- 
VB) channel model with the maximum delay spread 20 p,s 
and the number of paths p = 6 a was considered. 

Eig. d compares the mean square eiTor (MSE) performance 
of different channel estimation schemes. Both the static ITU- 
VB channel and time-varying ITU-VB channel with the mobile 
speed of 90 km/h in a 4 x 4 MIMO system were considered. 
The comb-type pilot based scheme used Np = 256 pilots, 
the TTOP scheme used Np = 64 pilots with T adjacent 
OEDM symbols for training where T = 4 for the time-varying 
channel and T = 8 for the static channel to achieve better 
performance, the TEJ scheme adopted time-domain training 
sequences of 256-length and Np = 64 pilots, and our proposed 
scheme used Np = 64 pilots with i? = 4 for fair comparison. 
Erom Eig.[3l we can observe that the conventional parametric 
TEJ scheme is inferior to other three schemes obviously. 
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Fig. 3. MSE performance comparison of different schemes in a 4 x 4 MIMO 
system, (a) Static channel; (b) Time-varying channel with the mobile speed 
of 90 km/h. 



Fig. 4. MSE performance of the proposed scheme in 4 x 4, 8 x 8 and 
12 X 12 MIMO systems. 


Meanwhile, for static ITU-VB channel, the MSE performance 
of the proposed parametric scheme is more than 2 dB and 5 
dB better than the TTOP and comb-type pilot based schemes, 
respectively. Moreover, for the time-varying ITU-VB channel, 
the superior performance of our proposed parametric scheme 
to conventional nonparametric schemes is more obvious. The 
existing sparse channel estimation scheme ||4| does not work, 
because path delays may not be located at the integer times of 
the sampling period for practical channels. The TTOP scheme 
works well over static channels, but it performs poorly over 
fast time-varying channels, since it assumes the channel is 
static during the adjacent OFDM symbols. Finally, the comb- 
type pilot based scheme performs worse than our proposed 
scheme, and it also suffers from much higher pilot overhead. 

Fig. in compares the MSE performance of the proposed 
scheme in 4 x 4, 8 x 8 and 12 x 12 MIMO systems. We can 
observe that the MSE performance of the proposed scheme 
in 12 X 12 MIMO system is superior to that in 8 x 8 MIMO 
system by 5 dB with the same Np, and outperforms that 
in 4x4 MIMO system with the reduced Np. These simulations 
indicate that with the increased number of antennas, the MSE 
performance improves with the same Np. Equivalently, to 
achieve the same channel estimation accuracy, the required 
number of pilots Np can be reduced. As a result, the total pilot 
overhead Vp_totai in our proposed scheme does not increase 
linearly with the number of transmit antennas Nt because 
the required Np reduces when Nt increases accordingly. The 
reason is that with the increased number of antennas, the 
dimension of the measurement matrix (e.g., H in ®) in the 


TLS-ESPRIT algorithm or the number of the sampling in the 
model of multiple filters bank Qol increases, thus the accuracy 
of the path delay estimate improves accordingly. 

The superior performance of the proposed scheme is con¬ 
tributed by following reasons. First, the spatially common 
sparse pattern shared among CIRs of different transmit-receive 
antenna pairs is exploited in the proposed scheme, such that 
we can employ the TLS-ESPRIT algorithm to obtain super¬ 
resolution estimations of path delays with arbitrary values. 
Meanwhile, the FRI theory indicates that the smallest required 
number of pilots is Np = 2P in a noiseless scenario. 
Therefore, the pilot overhead can be reduced as compared with 
conventional nonparametric schemes. Second, our scheme 
exploits the temporal correlation of wireless channels, namely, 
across several adjacent OFDM symbols, the sparse pattern of 
the CIR remains unchanged, and path gains are also correlated. 
Accordingly, by joint processing of signals of adjacent OFDM 
symbols based on (I8]l, the effective noise can be reduced 
and thus the accuracy of the channel estimation is improved 
further. 


V. Conclusions 

The proposed super-resolution sparse MIMO channel es¬ 
timation scheme exploits the sparsity as well as the spatial 
and temporal correlations of wireless MIMO channels. It 
can achieve super-resolution estimates of path delays with 
arbitrary values, and has higher channel estimation accuracy 
than conventional schemes. Under the framework of the FRI 
theory, the required number of pilots in the proposed scheme 
is obviously less than that in nonparametric channel estimation 
schemes. Moreover, simulations demonstrate that the average 
pilot overhead per transmit antenna will be interestingly re¬ 
duced with the increased number of antennas. 
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